Dynamics of recollisions for the double ionization of atoms in intense laser fields 
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We investigate the dynamics of electron-electron recollisions in the double ionization of atoms in 
strong laser fields. The statistics of recollisions can be reformulated in terms of an area-preserving 
map from the observation that the outer electron is driven by the laser field to kick the remaining 
core electrons periodically. The phase portraits of this map reveals the dynamics of these recollisions 
in terms of their probability and efficiency. 
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I. INTRODUCTION 

When subjected to short and intense laser pulses, the 
helium atom (or other atoms or molecules with two ac- 
tive electrons) may undergo double ionization The 
conventional route for double ionization is a sequential 
mechanism in which the field ionizes one electron after 
the other in an uncorrelated way. This process, called 
sequential double ionization (SDI) [lj, allows simple the- 
oretical predictions of double ionization yields : The dou- 
ble ionization probability is given by the product of the 
single ionization probability with the probability of ion- 
ization of the remaining ion. However, experiments car- 
ried out using intense linearly polarized laser fields 
have shown that in the range of intensity between 10 14 
and 10 15 W-cm~ 2 , double ionization yields depart signifi- 
cantly from the sequential predictions by several orders of 
magnitude. This observation has led to the identification 
of an alternative route to double ionization, called non- 
sequential double ionization (NSDI) [l| , in which the cor- 
relation between the two electrons cannot be neglected. 
Today NSDI is regarded as one of the most dramatic 
manifestations of electron-electron correlation in nature. 

Various mechanisms have been proposed to explain 
this surprise @, 0, [TTI j26|. When confronted with near- 
infrared experiments [aTl^. the recollision scenario [TP . 
[l2| seems in best accord with observations [2MT0I l2cl and 
is validated by quantum [l5l |29| s emi-classic al Il7l . l30l . 
HU and classical simulations [llM HI HI ■ Ac- 

cording to this scenario, a pre-ionized electron (referred 
as the "outer" electron [35[ ) , after picking up energy from 
the laser field, is hurled back at the parent ion by the laser 
and collides (thereby exchanging a significant amount of 
energy) with the remaining electron (referred as the "in- 
ner" electron (35[) trapped close to the nucleus. In gen- 
eral the inner electron experiences multiple recollisions 
and can eventually ionize, leading to double ionization 
if the outer electron remains ionized itself. Recollision 
has become "the keystone of strong-field physics" [l[ in 
the understanding of the electronic dynamics and light 
source design |36j. 

Even though the recollision mechanism is well settled 
in its broad outlines, some issues persist about the nature 



of collisions involved since not every recollision leads to 
double ionization (or does so right away). The efficiency 
of these collisions in transferring ionization energy during 
recollision have a direct bearing on the double ionization 
probabilities. 

In this manuscript, we investigate the dynamics of the 
recollisions which lead to double ionization, in particu- 
lar the energy exchange between the two electrons dur- 
ing successive recollisions. A number of NSDI features, 
obtained using classical models [ill E3, HH HI, l32| - [35| . 
are in very good agreement with results from quantum 
mechanical simulations and from experiments This 
agreement is ascribed to the prominent role of electron- 
electron correlation [l|, |2j [22| ■ In addition, classical me- 
chanical models have been used to a better understanding 
of the mechanisms because of their favorable scaling with 
system size. 

In what follows we consider the following Hamiltonian 
system describing, in the dipole approximation, a one- 
dimensional He atom using soft Coulomb potentials [2l|, 
|22| . l33l . |37| driven by a linearly polarized laser field of 
amplitude Eq and frequency lu : 



p 2 P 2 
H (x,y,p x ,p y ,t) = + 



1 



1 ' ^(x-yf + 1 
+ (x + y)E smut, 



(1) 



where x and y denote the position of each electron, and p x 
and p y their (canonically) conjugate momenta. The dura- 
tion of the pulse, i.e., the duration of the time-integration 
of trajectories, is 8 laser cycles. An analysis of typical tra- 
jectories of Hamiltonian (p} shows that the pre-ionized 
electron comes back to the inner region and exchanges 
energy with the inner electron several times which can be 
seen as repeated kicks delivered by the outer electron on 
the inner one (as seen in Fig.[T]). The key feature for dou- 
ble ionization is the energy exchanged during each kick. 
Viewing the recollision process as a periodic sequence 
of kicks suggests the use an area-preserving map which 
is constructed from periodically kicked dynamics, and 
widely used in various physical contexts in physics [38j . 
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The most prominent examples in atomic physics are the 
maps developed to model ionization of Rydberg atoms 
driven by microwave fields [H, HO] . 
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Laser cycles 



FIG. 1: Lower panel : Positions of the two electrons of 
a typical double ionizing trajectory of Hamiltonian {TJ as 
a function of time. The parameters of the laser field are 
I = 2 X 10 14 W ■ cm -2 and u = 0.0584 a.u. Upper panel : 
Interaction energy between the two electrons (green curve) 
defined as the soft Coulomb potential between the electrons. 
The dashed black line represents the threshold for recollision 
detection used in the analysis of recollisions. 

Discrete-time models for continuous-time periodic pro- 
cesses enjoy great popularity in physics. A paradigmatic 
example is the standard map [3a ] , the simplicity of which 
makes it an effective toy model for the study of chaotic 
properties in Hamiltonian systems. In addition, some 
rigorous properties can be derived, e.g., the transition 
from regular to chaotic behavior, the existence of ellip- 
tic periodic orbits, the persistence of rotational invariant 
tori, etc. Whether it is for integrating numerically such 
continuous-time processes or modeling physical phenom- 
ena, the main advantages of maps are that they can be in- 
tegrated more easily than continuous flows and that their 
properties appear with more clarity due to their reduced- 
dimensional phase space (as exemplified by Poincare sec- 
tions of continuous flows). Consequently they allow a 
deeper understanding of the dynamics and the underly- 
ing phenomena. 

Even though the nature of the kicks in the double ion- 
ization mechanism differs from what is observed for Ry- 
dberg atoms [39l . |4fJ , here we construct a map which is 
similar to the standard map. Using apt action-angle vari- 
ables for the inner electron dynamics, we construct the 
following map [4l| for the recollision dynamics : 

A n+ i = A n /(1 - esin^„), 

<Pn+i = (f n +Ty/2exp(aA n+1 ) + ecosifn, 

where A n and (p n are, respectively, the action and angle 
variables associated with the inner electron right before 



the n recollision. The constant a depends on the chosen 
atom, e.g., for He, a = — 9-\/2/16. The parameter e is the 
strength of the kick and will be related to the exchanges 
of energy at the recollision. The period of the kicks is 
denoted by T. The analysis of the phase portraits of this 
map shed some new light on the recollision dynamics. 

In Sec. HH we construct the map given by Eq. ([2]) from 
the analysis of recollisions experienced by the trajectories 
of Hamiltonian ([1]). In Sec. IIII| we analyze numerically 
this map in order to infer some properties of the recolli- 
sion dynamics and on the nonsequential route to double 
ionization. 

In what follows, instead of SDI and NSDI, we will use 
the more general UDI (uncorrelated double ionization) 
and CDI (correlated double ionization) (4l|. This ter- 
minology comes from the observation that a recollision 
may put the inner electron into an almost-bound state 
which then takes a significant time (sometimes more than 
one laser cycle) to ionize. With the previous definition, 
these so-called recollision excitation with subsequent ion- 
ization [42I l43l | events - by no means rare - would be 
labeled as SDI (because of the large time delay between 
the ionization of the two electrons) whereas they clearly 
correspond to a correlated process in the same way as 
NSDI. Here we consider CDI, where at least one recolli- 
sion is needed for double ionization, and UDI where no 
recollision is needed. 

Some of these results were anounced in a recent Let- 
ter [H|. 



II. DISCRETE-TIME MODEL FOR 
RECOLLISIONS 

Without the laser field (Eq = 0), typical trajectories 
associated with Hamiltonian ([IJ are composed of an elec- 
tron close to the nucleus (the "inner" electron) and one 
electron further away (the "outer" electron) [35|. This 
observation follows from the existence of four weakly hy- 
perbolic periodic orbits which organize the chaotic mo- 
tion. When the laser is turned on, the outer electron is 
quickly ionized while the inner one experiences a com- 
petition between the laser excitation and the Coulomb 
interaction with the nucleus. For nonsequential double 
ionizations, the inner electron is trapped in a bound re- 
gion about the nucleus and the only way to free itself 
is through a recollision with the outer electron when it 
returns to the core. We give an example of such a tra- 
jectory in Fig. Q] We note the fast first ionization of one 
electron while the other one remains trapped close to the 
nucleus. Because of the laser oscillations, the pre-ionized 
electron is hurled back at the core and recollides repeat- 
edly with the other one. After about 5 laser cycles, a 
final recollision manages to free both electrons and leads 
to a correlated double ionization. 

In a nutshell, the mechanism is the following one : 
When the outer electron returns to the core, the soft 
Coulomb interaction between the two electrons is no 
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longer negligible (particularly for the inner electron). It 
results in a rearrangement for this electron for which the 
action is modified. Since the outer electron comes back 
to the core quickly, this interaction is approximated by a 
kick experienced by the inner electron : When the outer 
electron returns to the core, it gives a kick in action to 
the inner one which jumps from one invariant torus to an- 
other (or to the unbound region, thereby ionizing). The 
action of the inner electron is constant between two rec- 
ollisions. 

In this section, we construct a simplified model for the 
recollisions which comes down to the following kicked 
Hamiltonian |4lj : 

N 

H m ((p,A,t) = H (A) + sAcos(f j^$(t-nT), (3) 

n— 1 

where Hq(A) is the integrable part of the Hamiltonian 
of the inner electron, and T is the delay between two 
recollisions, and e represents the strength of the kick and 
only depends on the intensity of the laser. Each part of 
Hamiltonian §3§ is designed from theoretical models and 
supported by statistical analysis of the recollisions. 

We construct an area-preserving map (From the kicked 
Hamiltonian given by Eq. ([3])). This is done in the 
standard way {38^ by defining ip n and A n as the an- 
gle and action of a trajectory of Hamiltonian Q at 
time (nT)~ (right before the n th kick). By integrating 
the trajectories between two kicks, i.e., from t = {nT)~ to 
t = [(n + 1)T] _ , we approximate the dynamics of Hamil- 
tonian ([3]) by the two-dimensional symplectic map ((2J. 

A. Integrable component : Ho {A) 

The effective Hamiltonian for the inner electron is 
given by (35[ : 

p 2 2 

#in (y,Py,t) = +yE smcjt, (4) 

1 vy + 1 

which is obtained from Hamiltonian (fTJ) by neglecting the 
interaction with the other electron. A quick inspection 
at its phase space shows that there are two main regions 
which result from the competition between the laser field 
and the Coulomb interaction : a bound region close to 
the nucleus where the electron remains bounded (on in- 
variant tori) at all times, and an unbound region where 
the electron leaves the nucleus quickly and ionizes [35j . 
The distance of the inner electron around the nucleus is 
best expressed in terms of its energy 

H (y,p y ) = ^-^=. (5) 

1 vy + 1 

The smaller this energy, the closer is the electron to the 
nucleus. In the absence of the field, Hq gives a natu- 
ral criterion for ionization, based on energy conservation 



for autonomous Hamiltonian systems. If the energy is 
smaller than zero, then the Coulomb interaction with the 
nucleus is strong enough to maintain the electron at a fi- 
nite distance for all times. On the contrary, if the energy 
is positive, then the electron will escape to infinity. 

In the neighborhood of the nucleus, the motion is har- 
monic with a frequency of y/2, and moving away from 
the nucleus, the frequency decreases. We observe numer- 
ically that the frequency v associated with H depends 
approximately linearly on the energy in the whole bound 
region (see Fig. [2J inner panel) : 

v (£) = a {£ + 2) + V2, 

where £ is the energy of the inner electron. We per- 
form a change of coordinates into action-angle in the 
system described by Hamiltonian ((SJ). From the equa- 
tion v(Hq) = 8H /dA, we obtain an expression for H : 

H (A) = -2 + V2?——±. (6) 
a 

The parameter a can be computed from a scries ex- 
pansion of the action associated with the inner electron 
around the bottom of the well. The equation for the 
action is 

where y m = y/A/£ 2 — 1 is the maximum position the 
electron can experience when it has energy £ . We define 
the energy as £ = —2^/1 — £ 2 , where £ is a small (posi- 
tive) parameter. Then, considering a series expansion in 
£, we end up with 

9^2 



Numerically, a ps —0.8 41]. In Fig. [5J wc compare the 
value for the action given by Eq. ^ with a numerical 
evaluation of the integral ([7]). 

B. The kicks 

1. Time delay between recollisions and number of 
recollisions 

The effective Hamiltonian for the outer electron is 

o 

P 

H out (x,p x ,t) = y + xE ° sincji ' 

which is obtained from Hamiltonian (fT]) by neglecting the 
interaction with the other electron and with the nucleus. 
Trajectories associated with Hamiltonian Hout are com- 
posed of a linear escape modulated by a sine function 
with the same period as the laser. It means that the 
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Action 



FIG. 2: Numerical evaluation of the energy given by Eq. ([5]) as 
a function of the action A (continuous red curves). The inset 
displays the frequency of the inner electron versus its energy. 
For comparison, in both panels, we give the values predicted 
by the approximate model © (dashed blue curves). 
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FIG. 3: Densities of recollision times (upper panel) and its 
Fourier decomposition (lower panel) obtained from a statisti- 
cal analysis of a large assembly of trajectories for an intensity 
I = 2 x 10 14 W ■ cm" 2 and laser frequency cj = 0.0584 a.u.. 



outer electron typically experiences two returns to the 
core per laser cycle. 

We have collected statistical data from 16000 trajecto- 
ries associated with Hamiltonian ([TJ at recollision times, 
starting with a microcanonical initial distribution over 
the ground state energy surface H(t = 0) = —2.24 
a.u. [211 . [35l Numerically looking at the peaks in 

the energy of interaction between the electrons (defined 
as the soft Coulomb potential between the electrons) re- 
veals times of recollision (see upper panel of Fig.Q]). From 
there on, one can collect and analyze characteristic data 
of the recollision process such as times when they take 
place, momentum of the outer electron, number of rec- 
ollisions and exchanged energy (or action) during non- 
ionizing recollisions. One can define an energy for the 
inner electron, defined by Hq, as long as it is trapped in 
the bound region; it results in the impossibility to mea- 
sure the amount of exchanged energy during recollisions 
leading to ionization. Consequently, ionizing recollisions 
are systematically discarded from the statistical analy- 
sis. In Fig.[3J we give an example of the density of return 
times of the outer electron and its associated spectral 
decomposition. It reveals that the main frequency for 
recollisions peaks around 2 per laser cycle which corre- 
sponds to a main period between two recollisions of half 
a laser cycle. As a result, in the map we set the time 
delay between successive kicks to be equal to T = n/u>. 
Since the time duration of the laser pulse is 8 laser cycles, 
the inner electron experiences at most 15 recollisions. In 
Fig. |H we display the statistical distribution of the num- 
ber of recollisions collected from the analysis of typical 
trajectories and typical double ionizing ones, for a fixed 
intensity. It shows that most of the trajectories do not 



undergo any recollision and typically between 2 to 4 rec- 
ollisions are required to trigger correlated double ioniza- 
tion. The density depends weakly on the intensity in the 
intermediate range of intensity. It indicates that the map 
should not be iterated more than 12 times to reproduce 
accurately the dynamics of the recollisions experienced 
by the inner electron. 




Number of recollisions 



FIG. 4: Distribution of the numbers of recollisions for a laser 
intensity I = 2 X 10 14 W • cm" 2 and frequency u = 0.0584 
a.u.. Dark blue bars correspond to data collected from typical 
trajectories (including non-ionizing ones) and light green bars 
correspond to data associated with double ionizing trajecto- 
ries. For a better layout, the middle part of the statistics has 
been cut for the zero recollision component which goes up to 
70%. 
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2. Exchanged action during recollisions 

In Hamiltonian ([3]), the recollisions are modeled by a 
kick in action equal to eAcosy such that a kick might 
increase or decrease the action according to the respective 
phase between the two electrons. In addition, it is more 
difficult to kick the inner electron out if it is at the bottom 
of the well, so the kick strength is proportional to A. In 
this way, the action remains positive at all times since 
A = is invariant. The maximum strength of the kick 
depends strongly on Eq. 

It is well-known that the maximum energy the outer 
electron can bring back to the core is equal to £ ma x = 
tdJ p where U p = E%J (4a; j is the ponderomotive energy 
and k ~ 3.17 [TJ [12, Eljj- However, an inspection of 
the energy exchanged during recollisions shows that this 
amount is significantly smaller (see Fig. [5]). Below, we 
analyze the recollisions in order to explain the trends 
observed in Fig. [5l 




/ (Wcm ) 



FIG. 5: Standard deviation (red stars, left hand vertical scale) 
of shared action during recollisions as a function of the laser 
intensity /. An approximation of the standard deviation (red 
continuous curve) is given by Eq. (|l(jp . The thin black contin- 
uous line corresponds to the maximum recollision energy £ max 
(expressed in terms of actions) and the vertical dotted line in- 
dicates the intensity after which we stop iterating the map. 
The dotted curve (blue line, right-hand vertical scale) shows 
the mean number of recollisions computed from all analyzed 
trajectories. 

a. Low intensity limit For an inner electron at the 
bottom of the well (A — 0) which experiences a kick 
with energy kU p , its net change in action is thus equal to 
AA = log (l + KUpa/V2) /a, using Eq. ©. To leading 
order it gives 



AA: 



for small Eq. It confirms the increase as Eq observed 
through the analysis of recollisions when the amplitude 



of the laser is weak (see Fig. [5]). 

b. High intensity limit As the intensity of the laser 
field increases up to 10 14 W • cm -2 , the energy the outer 
electron brings back to the core (and potentially to the in- 
ner electron) varies approximately as Eq. Increasing the 
intensity further, it appears from the analysis of recolli- 
sions that the exchanged action decreases after a critical 
intensity of 2 x 10 14 W ■ cm~ 2 . To model the process of 
recollision, we consider a Hamiltonian for the inner elec- 
tron (y) when the outer one (ar) comes back to the core 
as : 

H YCC ( 



(y,p y ,t) = — 
1 

^{y-x{t)f^ 



l 



= + yEo sin tut, 



(8) 



1 



which is obtained from the reduced Hamiltonian ((4]) of 
the inner electron by adding a passive soft Coulomb inter- 
action with a quickly returning outer electron. Because 
of its large momentum, we consider that the outer elec- 
tron is not affected by the interaction with the inner one 
(i.e. we impose the dynamics for the outer electron in- 
dependently of the dynamics of the inner one) and the 
nucleus. From its position and momentum at recolli- 
sion, respectively yo and po, its trajectory is almost a 
straight line : x (t) = yo + p^t, where we have set the 
origin of times at the time of the recollision. For large 
intensities, the inner electron in the bound region is close 
to the nucleus, since the bound region becomes smaller, 
and in such a configuration the dynamics is well approx- 
imated by harmonic potentials. Finally, we assume the 
motion for both the inner and outer electrons to be much 
faster than the laser field variation. Thus it is assumed 
that locally Egsinuit « tjEq. As a result, Hamiltonian 
.ffrec given by Eq. ([5]) is simplified through an expansion 
around the position y = as 



•ffrec {y,Py,t) 



P 2 y 



y 



r)E y- (9) 



2 " 2 

This model of recollision is valid as long as the two elec- 
trons are close enough to each other to allow an effec- 
tive exchange of energy. We denote by L the maxi- 
mum distance between them to have an effective inter- 
action. When the outer electron is further away, the 
term (y — x) 2 is canceled out. The resulting interval of 
time (— r, t) during which the two electrons are interact- 
ing with each other is equal to r = L/po. When the 
outer electron leaves the region of interaction, the effec- 
tive model for the inner one becomes : 



Pi 



y 



v E oy- 



The resulting exchanged 
lision is equal to AH 



energy during the recol- 
= H har (y(T),p y (T)) - 



-ffhar {y (— t) ,p y (— t)). The positions y (±t) and mo- 
menta p y (±r) are computed by (forward and back- 
ward) integration of the trajectory with initial condi- 
tions y (t = 0) = yo and p y (t — 0) = po (the inner 
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electron momentum at recollision) and whose dynam- 
ics is given by Hamiltonian ©. We expand these ex- 
pressions up to order 0(r 4 ), and the leading term of 
AH is g iven by 2pqT)Eqt 3 /3. From the recollision pic- 
ture [ill. l45j the maximum momentum the outer electron 
can have when it returns to the core is pg nax = K 'Eo, 
where n' = \J k/(2uj 2 ). Statistical analysis of trajectories 
shows that for large intensities, the momentum distri- 
bution of the outer electron is centered around ±7'po lax 
where 7' « 0.8 (see Fig. [5]). As a result, we choose the 
outer electron's momentum proportional to the field am- 
plitude : po — jEq (i.e. 7 = j'k') and the trajectory of 
the inner electron during recollision can be computed an- 
alytically. It allows one to consider a series expansion for 
the exchanged energy for the inner electron. From AH, 
we recover the net variation in action at the recollision 
through Eq. @ : 



AA = 



AH 

7f 



OH 2 



To leading order, the net exchange in action is equal to 



AA 



V2i]L 3 1 
3 7 2 £0 



O 



1 

^0 



This simple model of recollision explains the decrease as 
1/Eq for the action exchange observed during the anal- 
ysis of recollision, provided that the interaction length 
L is independent of Eq. Moreover, A A varies as u 2 , 
which means that we expect less exchanged energy at 
low frequency for high intensities. So the correlated dou- 
ble ionization probabilities are expected to be lower in 
the high intensity regime, and higher for low intensities 
as the laser frequency decreases. 

In summary, to combine the two trends of the mean 
shared action AA (Eq) (proportional to Eq at low inten- 
sities as given by and to I/Eq at higher ones) we fit 
it by : 



AA (E 



olE 2 
1 + ml ■ 



(10) 



The parameter a is equal to k/(4v / 2w 2 ). The parameter 
P is given by /3 = 3k 2 7 /2 / (16t]L 4 (jj 2 ) and is obtained by 
a numerical fit so as to accurately reproduce the evolu- 
tion of the mean exchanged action during recollisions (see 
continuous lines in Fig. [SJ. For instance, for oj = 0.0584 
a.u., the fitted value for /3 is = 5.1 x 10 3 . We consider 
for the kick strength in our map, e = A A (Eq) as given 
by Eq. (JTOJ). 



III. NUMERICAL ANALYSIS OF THE MAP (HJ) 
A. Inner electron distribution 




FIG. 6: Outer electron momentum distribution at recollision 
times as a function of the intensity of the field. The density is 
obtained from the analysis of a large assembly of trajectories 
of Hamiltonian {TJ for cj = 0.0584 a.u.. For comparison we 
also display momenta po — ij'n'Eo where 7' = 0.8 (blue 
line). 



an outer electron : the inner one is the electron with the 
smallest energy. We compute the distribution in action 
for the inner electron as defined in Sec. Ill Al An inspec- 
tion of the shape of this density reveals that it has an 
exponential decrease on the range of accessible actions 
(see Fig. [7]). To model this distribution, we choose an 
exponential law with a truncated tail given by 



f,. An {A) 



A 



1 



3 -AA 



(11) 



where Aq is the maximum allowed action for the inner 
electron and Xn ( x ) = 1 if 2; € and otherwise. The 
two parameters A and Aq are adapted so as to agree 
with the distribution deduced from the microcanonical 
set. The first one, Aq, is the maximum allowed action 
for the inner electron and it is computed from Hamilto- 
nian ([1]) for Eq = 0. For the parameter A, we consider 
a numerical fitting with a large assembly of initial con- 
ditions on the ground state energy surface. The fitting 
is done using the maximum likelihood method [46j (the 
mean value method gives the same result). The numer- 
ical evaluation of the two parameters yields Aq = 0.61 
and A = 5.3. In Fig. [7] we compare densities in action for 
the inner electron obtained from a microcanonical initial 
distribution and the exponential law given by Eq. (jlll) . 
Note the very good agreement between the two distribu- 
tions for almost all allowed actions. 



B. Phase portrait of map ([2]) 



Starting from a microcanonical distribution on the 
ground state energy surface one can identify an inner and 



Now that the parameters of the map as well as the 
initial conditions are determined, we investigate numeri- 
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FIG. 7: Inner electron distribution function versus action. We 
compare the density obtained from a microcanonical set of 10 5 
initial conditions (large black curve) with the exponential law 
(thin red curve). For comparison, we display the initial 
maximum action for the inner electron Aq (dashed red line). 



cally the dynamics given by map @. In Fig. [51 we display 
two phase portraits for two laser intensities : One in the 
intermediate range of intensity (I = 2 x 10 14 W • cm" 2 ) 
shows a phase portrait which appears to be very chaotic, 
and one at high intensity (I — 2 x 10 15 W • cm" 2 ) shows 
a more regular phase portrait. In the chaotic region, the 
diffusion is much stronger in the intermediate range of in- 
tensities than for larger intensities (see Fig. |5J left panel, 
where trajectories quickly escape from the core region, 
explaining why there are fewer points than in the right 
panel) . Since the strength of the kicks decreases with the 
intensity at high intensity, the phase space becomes more 
regular. If the inner electron is inside an elliptic island 
(which occurs mainly at high intensities), it will not ion- 
ize regardless of the number of recollisions it undergoes. 
As the intensity increases, the recollisions become less 
effective and the map becomes integrable so fewer CDI 
events occur. Therefore there are two competing mech- 
anisms for the vanishing of the CDI probability at high 
intensity : the decrease of the size of the bound region, 
and the lack of efficiency of the recollisions due to the 
regularity of the dynamics of map ([2]) . 



C. Statistical analysis versus intensity 

Through map |(3J), we have derived a simple model 
for the dynamics of recollisions experienced by the in- 
ner electron initially in the bound region. From this 
model, we compute ionization probabilities of the in- 
ner electron from which we deduce the double ionization 
probability : Double ionization occurs when the inner 
electron ionizes provided we assume that the outer elec- 



FIG. 8: Phase space portraits of some trajectories of the 
map @ with ui = 0.0584 a.u. for low intensity (left panel) 
I = 2 x 10 14 W • cm" 2 , and for high intensity (right panel) 
I = 2 x 10 15 W • cm -2 , represented by vertical dashed lines 
in Fig. In the right panel, we indicate the critical ac- 
tion A m = 0.71, after which the inner electron ionizes, by 
a horizontal line (whereas the critical action A m is 1.60 in the 
left panel). 



tron remains ionized for all times. The picture of the 
bound and unbound regions for the effective Hamilto- 
nian ((4]) of the inner electron gives a natural criterion 
for ionization : Once the inner electron has reached an 
action larger than the outermost invariant torus (with 
action A m ), it is driven away from the nucleus by the 
laser field. Therefore, all recollisions leading to an ac- 
tion larger than a critical value A m (which depends on 
Eq) subsequently lead to ionization of the inner elec- 
tron. In angle-action variables, the unbound region be- 
comes T>(E ) = {(ip,A) s.t. A> A m (E )}. We repre- 
sent A m (Eq) as a function of the laser intensity in Fig.HUl 
From an initial distribution in angle-action coordinates 
obtained from the microcanonical distribution (jlll) of in- 
ner electrons in phase space, we iterate map @ a fixed 
number of times for different intensities (and thus differ- 
ent e). In what follows, we compare the results given by 
the map to the probability obtained using a direct inte- 
gration of the trajectories of Hamiltonian (JTJ. As a func- 
tion of the intensity of the field, these probability curves 
take the form of a "knee" 0-0, EM, El HI HI El which 
shows an enhancement of the double ionization probabil- 
ity in the intermediate range of intensities. In Fig. |H1 we 
display the double ionization probability as a function of 
the laser intensity. We disregard recollisions for intensi- 
ties larger than 2.5 x 10 15 W • cm" 2 . This adjustment 
is motivated by the weak probability of recollisions we 
have detected in the data analysis (see Fig. [5]). A quick 
inspection of Fig. [9] reveals a bell-shaped curve for the 
resulting nonsequential component (35l . |4lj . We notice 
that it qualitatively reproduces the trends observed in 
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the double ionization yields observed using a statistical 
analysis of trajectories of Hamiltonian ([T]) for two differ- 
ent values of the laser frequency. In particular, the asym- 
metry in the increase and decrease of the nonsequential 
component is worth noting. 

A rather intuitive mechanism to explain the decreas- 
ing part of this bell shape is the conversion from nonse- 
quential trajectories into sequential ones when the laser 
field becomes stronger. However, in Fig. [9l we notice 
a local decrease of the total yield (also observed with 
quantal computations [13]) which is larger by several or- 
ders than the increase of sequential process. We notice 
that this incompatibility is readily observed in Fig. 1 of 
Ref. [13] . As explained previously (see Sec. Ill B 2b|) . the 
decrease of the nonsequential component is mainly due to 
the decrease of recollision efficiency with the laser inten- 
sity. At large intensity, the UDI probability is given by 
the proportion of the ground state energy surface where 
both electrons belong to T>(E ). In Fig. |H1 we display 
the UDI component predicted by this model (continuous 
blue curve) which is in good agreement with the dou- 
ble ionization probability obtained by integrating the full 
Hamiltonian ([1} in the high-intensity regime. 



IV. CONCLUSION 

In summary, we have studied key properties of the in- 
elastic electron-electron recollisions through the analysis 
of recolliding trajectories. We connected our findings to 
a simplified model for the dynamics of the recollisions, 
which amounts to an area preserving map in the action- 
angle coordinates of the inner electron. The statistical 
analysis of the (discrete-time) trajectories of this model 
results in the hallmark "knee" shape for the probability 
of double ionization versus intensity. A proper decompo- 
sition into correlated and uncorrelated double ionization 
yields a bellshape for the correlated process and a mono- 
tonic rise for the uncorrelated one. 
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FIG. 9: Double ionization probability for Hamiltonian |T]) (continuous lines) and the map ([2| (red dark circles) as a function 
of the laser intensity / for two laser frequencies (left panel uj — 0.0584 a.u. and right panel oj = 0.02 a.u.). Because of the low 
probability of recollision, we stop iterating the mapping after / = 2.5 x 10 15 W-cm~ 2 for lj — 0.0584 a.u. and I = 10 15 W ■cm 
for ui = 0.02 a.u. For each frequency, we also indicate the expected UDI probability obtained as described in the text (green 
light circles) . Vertical dashed curves refer to intensities where phase portraits are displayed in Fig. [5] 
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